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The paper presents performance and robustness analysis of the modified reference model 
MRAC (model reference adaptive control) or M-MRAC in short, which differs from the 
conventional MRAC systems by feeding back the tracking error to the reference model. 

The tracking error feedback gain in concert with the adaptation rate provides an additional 
capability to regulate not only the transient performance of the tracking error, but also 
the transient performance of the control signal. This differs from the conventional MRAC 
systems, in which we have only the adaptation rate as a tool to regulate just the transient 
performance of the tracking error. It is shown that the selection of the feedback gain and the 
adaptation rate resolves the tradeoff between the robustness and performance in the sense 
that the increase in the feedback gain improves the behavior of the adaptive control signal, 
hence improves the systems robustness to time delays (or unmodeled dynamics), while 
increasing the adaptation rate improves the tracking performance or systems robustness 
to parametric uncertainties and external disturbances. 

I. Introduction 

Robustness analysis of classical linear time invariant (LTI) systems are based on the notions of gain and 
phase margins, which provide stability guarantee in the presence of unmodeled dynamics and unstructured 
uncertainties. These stability metrics are computed using linear analysis tools, which cannot be easily 
extended to inherently nonlinear adaptive system. While the notion of the gain margin for the adaptive 
systems can be essentially defined in the same way as for the LTI systems, and the estimate of the gain 
margin can be obtained using the Lyapunov analysis tools (see for example [12]), the notion of phase margin 
is not defined for the nonlinear systems. On the other hand, time delay margin has been accepted as a viable 
alternative for the robustness analysis of adaptive systems. 

When the adaptation mechanism is turned off, the controller can be viewed as a fixed gain controller, 
and the adaptive system reduces to an LTI system, when the underlying dynamics is linear. Therefore, it is 
tempting to freeze the adaptive parameters over a time window, in order to apply the linear analysis tools 
to estimate the stability margins of the system. However, there is no guarantee that the obtained margins 
in this fashion relate to the stability of the actual adaptive system. Recently, some methods have been 
proposed for approximate phase and gain margin analysis of the adaptive systems without terminating the 
adaptation. In [11], it is proposed to define an LTI system bounding the closed- loop adaptive system then 
locally evaluate the phase and gain margins for the bounding LTI system. In [14], an LMI based approach 
is proposed to calculate the stability margins by transforming the nonlinear adaptive system into a linear 
parameter varying system. 

In [9], several linear adaptive control systems in the presence of measurement noise and external distur- 
bances have been compared using the linear analysis tools from the point of view of the performance and 
robustness. 

The time delay margin of the adaptive systems are not readily computable as well. Some methods have 
been reported in the literature which use Lyapunov based stability tools to estimate it. In [1], the authors 
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use Pade approximation to transform the original time delay system to a higher order system without time 
delay and approximate the time delay margin. In [7], the Razumikhin method has been use to iteratively 
estimate the time delay margin. These methods are general to any adaptive control systems, but generate 
conservative results. In [3] a lower bound for the time delay margin has been obtained that is applicable 
only to the C\ adaptive control system. 

In this paper we do not pursue a goal to estimate the time delay margin of the nonlinear M-MRAC 
systems, but to show how can one improve it by the proper choice of the design parameters. For this 
purpose, we follow the Razumikin’s approach (see for example [5, 10]) to show that the robustness of M- 
MRAC to time delays can be improved by increasing the error feedback gain, hence the tradeoff between the 
performance and robustness can be resolved by operating the adaptation rate and the error feedback gain. 
It is also shown that the M-MRAC design has improved robustness to unmodeled dynamics compared to 
conventional MR AC design. 

In addition to time domain analysis, frequency domain analysis is conducted for the linear M-MRAC 
systems. Using methods from linear systems theory, several critical transfer functions are computed between 
the systems input and output signals. These transfer functions, which are called Gang of Six in [2], allow us 
to compute the stability margins and analyze the performance of the close-loop system in the presence of the 
measurement noise and external disturbances. By explicitly computing the crossover frequencies and gain 
and phase margins, it is shown that: 1) the gain margin increases when the error feedback gain is increased 
and decreases when the adaptation rate is increased; 2) the gain crossover frequency decreases when the 
feedback gain is increased and increases when the adaptation rate is increased, 3) the phase margin and 
the time delay margin increase when the feedback gain is increased and decrease when the adaptation rate 
is increased. On the other hand, the tracking performance and disturbance attenuation properties of the 
closed loop signal are improved with the increase in the adaptation rate. Hence, the same tradeoff between 
the performance and robustness is obtained as in the time domain analysis. 

The rest of the paper is organized as follows. In Section II we give preliminaries on the M-MRAC design. 
We provide the robustness analysis to time delays in Section III. The frequency domain analysis of linear 
M-MRAC systems is given in Section IV. In Section V the robustness of M-MRAC to unmodeled dynamics 
is provided. Some concluding remarks are given in Section VI. 

II. M-MRAC 

Consider a multi-input multi-output uncertain linear system 

x(t) = Ax(t ) + Bu(t), :r(0) = x G , (1) 

where x £ K™ and u £ R 9 are the state and input of the system respectively, and A £ R™ xn and B £ R nX9 
are unknown constant matrices. The control objective is to track a reference model 

x r {t~) — A m :r r (t) -f- , 3? r (0) = Xq , (2) 

where A m £ R nx ™ is a Hurwitz matrix and B m £ R nx? is of full column rank, assuming that there exists a 

matrix K\ £ R” x<? and a sign definite matrix A £ R ?x? such that the matching conditions hold 

B = B m A (3) 

A = A m — BKj . 

The external command r(t) is usually assumed to be bounded and piece-wise continuous, which has a 
bounded derivative almost everywhere on [0, oo). Therefore, it is piecewise uniformly continuous on [0, oo) 
in the sense of the definition from [15]. We notice that discontinuous commands such as steps or square 
waves satisfy these conditions. 

The system (1) can be written in the form 

x(t) = A m x(t) + B m r(t) + B m A[u(t) - 0 T /(t)] , (4) 

where 0 = [Kj Kj], f(t ) = [x T (t) r T (t)] T and Kj = A” 1 . This suggests that the control signal can be 
defined as 


u{t) = Q r (t)f(t) 


( 5 ) 



where 0(f) is the estimate of the unknown parameter 0, and is updated online according to the adaptive 
law 

©(*) = - i yf{t)e? n (t)PB m , 0(0) = 0 O , (6) 

where 7 > 0 is the adaptation rate, and P = P T > 0 is the solution of the Lyapunov equation 

A T m P + PA m = -Q (7) 

for some Q = Q T > 0. The error signal e m (t) = x(t) — x m (t) is defined with respect to modified reference 
model 

— A m x m (t) + B. m v(t) T Xe m (t), 3? m (0) xq , (8) 

where A > 0 is a design parameter. 

Introducing the parameter estimation error as 0(f) = 0(f) — 0, the following error dynamics can be 
derived 

e m (t) = (Am - AI)e m (f) + B m AQ T (t)f(t) . (9) 

It is easy to see that the actual tracking error eft) = x(t) — x r (t) satisfies the equation 

^[e m (f) - e(t)] = A m [e m (t ) - e(t)] - A e m (t) , (10) 

since e m (t) — eft) = x r (t) — x m (t). We notice, that setting A = 0 results in conventional MRAC design with 
^m(^) — ^(t)- 

The following theorems give the asymptotic and transient properties of the M-MRAC architecture in the 
ideal case. 

Theorem 1 Let the system (1) be controlled by the M-MRAC scheme given by (5), (8) and (6). Then all 
closed-loop signals are bounded, and e(t) —> 0 , e m (t) —t 0 and u(t) — > 0 as t — > 00 , where u(t) = Q T (t)f(t). 

The theorem is proved using the Lyapunov function 

V(t) = em(t)Pe m (t) + ^ trace ^A0 T (t)0(t)^ (11) 

and the extended Barabalat’s lemma form [15]. 

Theorem 2 Let the system (1) be controlled by the M-MRAC scheme given by (5), (8) and (6) and 

A = 2a'yX max (B^ l PB m ) . (12) 

Then the following bounds hold 

IKBII < fre"* + &7- 1 (13) 

l|e(i)|| < 

where a = || a^o II 2 + sup t ||r(t)|| 2 , fti,f32,P3 and /I4 are positive constants independent of 7, v is proportional 
to ypy, v* = nrin[i/, v m \ and v m is the rate of decay of e Amt . 

The proof of the theorem can be found in [13]. 

III. Robustness to Time Delays 

In this section we analyze the properties of the M-MRAC design in the presence of a time delay in 
the input signal. In this case, the system is given by the so called retarded functional differential equation 
(RFDE), for the analysis of which we use Razumikhin’s method. Here, we give some background information 
on the method. The interested reader can get more detailed information in the monographs [5] and [10]. 



A. Preliminaries on Razumikhin’s method 

Let the system be given by a retarded functional differential equation (RFDE) (the terminology and defini- 
tions are adopted from [5] and [10]) 


z(t) = , (14) 

where z t (6) = z(t + 6) for —h < 9 < 0 for a given real number h > 0, and z to = <p. Here, / : I x C[—h, 0] — > 
R” is a continuous functional which is locally Lipschitz in the second argument, <f> £ C[—h, 0], and C[—h, 0] 
denotes a Banach space of ?r-dimensional real- valued vector-functions with the supremum norm. 


Definition 1 ( [5], p. ISO) Let f(t, z t ) = 0 for all t £ [to, oo). The solution z = 0o/ the RFDE (If) is said 
to be uniformly stable if for any to and any e > 0 there exists a 6 = 6(e) > 0 such that <j> £ B(0,<5) implies 
that Zt(to,4>) £ B(0, e) for all t, > to, where z t (to,4>) solution of the RFDE (If) through (to,<j>), and 

B(0,r) denotes a ball of radius r centered at zero. 

One way to establish the uniform stability of the trivial solution of RFDE (14) is to check the conditions of 
the Rasunikhin’s theorem (see [5], p. 152) 

Theorem 3 Suppose u,v,w : R + —> M + are continuous, nondecreasing functions, u(s),v(s) positive for 
s > 0, u(0) = u(0) = 0, v strictly increasing. If there is a continuous function V : R x R” — > R such that 


and 


If 


«(NI) < V(t,z) < v(||z||), 
V(t,z(t)) < — w(||z(t)||) , 
V(t + 9, z(t + 9)) < V (t, z(t )) 


(15) 

(16) 
(17) 


for 6 £ [— h,0], then the solution z = 0 of the RFDE (If) is uniformly stable. 

As it is shown in [10], the stability of the RFDE (14) can be analyzed using the same Lyapunov function 
used for the analysis of the ordinary differential equation resulting from RFDE (14) when h — 0. 


B. Nonlinear analysis 

To analyze the robustness properties of M-MRAC system we insert a time delay td at the input point of the 
system. That is the control signal is defied as 


which results in the system 


u(t.) 


0, 0 <t<t. d 

u(t-t d ), t>t d 


Obviously, we have 


x(t) 


Ax(t), 0 < t < t d 

Ax(t) + Bu(t - td), t>t d 


x(t) = e At x 0 


(18) 


(19) 


( 20 ) 


on the interval 0 < t < t d . The corresponding modified error dynamics take the form 


L-m(t) 


(A m - AI)e m (f) + (A - A m )x(t ) - B m r(t), 0 < t < t d 

(A m - AI)e ro (t) + B m A<d T (t)f(t) - B m A u(r)dT, t > t. d 


( 21 ) 



The error signal e m (t) is readily computed on the delay interval 0 < t < t d to be 

e m (t) = (A- A m ) ( e (^m-Ai)(t-r) | 'e Ar Xo + B m r(r)] dr . (22) 

Jo 

The stability of the error system (21) along with the 0(f) equation readily, which is obtained from the 
adaptive law (6), is analyzed using the Lyapunov function V (11). First of all, we notice that 

A mi „(P)||e m || 2 + ^^||0|| 2 F < V(e m , 0) < A max (P)||e m || 2 + Amax(A) ||e||| , (23) 

7 7 

where the subscript F denotes the Frobenius norm. Obviously the inequality (15) holds with u(z) = 
Amin (P) || e m || 2 + Ami ° (A) ||0 |If and v{z) = A max (P)||e m || 2 + Amax(A) ||0|| F . Next, we evaluate V for t > t d . It 
is straightforward to check that 


V(t) = -e^(t)Qe m (f) - 2A e^(t)Pe m (t) - 2e^(t)PP m A f u(t)cLt . (24) 

Jt-t d 

Our goal is to investigate the negativity of V (t) depending on the design parameters A and 7. To this end, 
we compute u(t) using the adaptive laws and the system’s dynamics. 

u{t) = —'yB^ n Pe m (t)[x 2 (t) + r 2 (t)\ + [Kj A + BkJ]x(t) + kj BKjr(t) + kjr(t) . (25) 

Therefore V(t) (25) can be expressed as 

V(t) = — (t)Qe m (t) - 2Ae^(t)Pe m (t) + 2jai(t,td) + s 2 (t,t d ) . (26) 

where we denote p{r) = x 2 (t ) + r 2 (r), 


si{t,td) = 2e^(t)PP TO APj)P f e m (r)p(r)dr 

Jt-t d 

si(t,t d ) = 2e^(f)PP m A f {[Kj ( t)A + BKj (r)]£c(r) + kj ( t)BKJ ( r)r(r ) + kj (r)r(r)}dr 
Jt-t d 

The conventional MR AC design corresponds to A = 0. In this case the design parameters are Q and 
7, assuming that A rn and B m are given from the performance specifications. The matrix Q can be chosen 
from the point of view of maximizing the ratio A m i n (Q)/A max (P), which results in Q = I (see [8]), p. 342 
for details). Hence, 7 is the only design parameter to be selected from the perspective of minimizing the 
Lyapunov derivative 

V(t) = -e^(t)e m (t) + 27Si(f,fd) + s 2 (t,t d ) , (27) 

which can be upper bounded as follows 

V(t) < -||e m (t)|| 2 + 27t d ||PP m AP^P||||e t (r)||| oo ||p t (r)|| z ; 00 + s 2 (t,t d ) . (28) 

Here, ||e t (r)||£ ao is the norm in the extended space [8] (p. 197). Clearly, the upper bound on V(t) 
increases with respect to 7. In fact, it is increasing in the product 7 t d , which implies that the time delay 
margin is adversely affected by the large adaptation rate. This is consistent with practical applications and 
linear analysis (see for example [4] and references therein). 

On the other hand, for the M-MRAC design the upper bound (28) has the form 

V(t) < -||e m (t)|| 2 - 2AA min (P)||e m (t)|| 2 + 27t d ||PR m AB^P||||e t (T)||£ oo ||p t (r)||£ 00 + s 2 {t,t d ) , (29) 

which implies that the upper bound on V(t) is a decreasing function in A. This observation gives as a 
capability to increase the time delay margin of M-MRAC systems by the proper choice of A. It follows from 
the Razumikhin’s theorem that if for given initial conditions and design parameters Q = I, 7 and Ai the 



maximum time delay, for which the error system ( 21 ) is stable, is t d i , then for some A 2 > Ai the inequality 
td 2 > tdi holds, where the error system ( 21 ) is stable for the time delay t d 2 - 

We test this assertion in a simulation example for a first order system with A = 2, B = 1 , A m = — 1 and 
B m = 1. We choose a unit step at time t = lsec. as an external command, set 7 = 1000 and compute A 
according to equation (12). Figure 1 displays the performance of conventional MRAC (A = 0) and M-MRAC 
for xo = 0. It can be noticed that the MRAC performance is not satisfactory even in the ideal case. Next 
we insert a time delay td at the input point and numerically compute the system’s time delay margin. The 
conventional MRAC goes unstable at td = 0.016sec. The time delay margin of the M-MRAC is presented in 
Table 1, where A* = 2^/y, since xq = 0 and supr 2 (f) = 1 . 


Output time history 



Input time history 





Output time history 





(a) MRAC 


(b) M-MRAC 


Figure 1. Input and output time histories in the ideal case. 


Table 1. Time delay margin variation with respect to A. 


Error feedback gain A 

A* 

2A* 

4A* 

Time delay margin 

0.039 

0.80 

0.126 


Obviously the system’s robustness to time delays is increased with the increase of A, which confirms the 
analytical predictions. The system’s performances corresponding to three cases in Table 1 are displayed in 
Figure 2. It can be observed that while the control objective is met and the system’s output converges to 
reference model’s output for all values of A, the transient performance degrades when A is increased from the 
suboptimal value A* . That is A can be used to resolve the tradeoff between the robustness and performance 
for the given adaptation rate 7 . 



(a) A = A*, t d = 0.039 



(b) A = 2A*, t d = 0.080 



(c) A = 4A*, t d = 0.126 


Figure 2. System’s input and output time histories with variable feedback gain and time delays. 


It has been noticed that setting the error feedback gain to a time variant quantity A* ( t ) = 2-^/y [x 2 (t) + r 2 (t) ] 
improves the both the performance and robustness to time delays of the M-MR AC system. In this case The- 


orem 1 is still valid, hence the above reasonings based on the Razumikin’s theorem apply. Table 2 presents 
the time delay margins of the system with variable A *(f), which are substantially improved over the fixed A 
case presented in Table 1. 


Table 2. Time delay margin variation with respect to time varying A. 


Variable gain A (t) 

A *{t) 

2 A*(f) 

4A*(t) 

Time delay margin 

0.055 

0.104 

0.162 


The corresponding input and output time histories are displayed in Figure 3 



(a) A = A „(t), t d = 0.055 



(b) A = 2A*(i), t d = 0.104 



(c) A = 4A*(t), t d = 0.162 


Figure 3. System’s input and output time histories with variable feedback gain and time delays. 

As in the case of fixed A, control objective is met in all cases, but the system’s performance degrades 
with the increase of A at the expense of improved robustness to time delays. 


C. Linear analysis 

When Ki = 0 and r = const the M-MRAC design results in the linear time invariant system. In this case, 
the adaptive control is just an integral control. The system can be written in frequency domain as 

x(s) = (si — A m )~ 1 Bu(s) 

x m (s) = (si - A m ) _ 1 Rr(s) + A(sl- A m ) _ 1 [a;(s) - cc m (s)] 

u(s) = -B^P-[x{s) - x m (s)\, (30) 

s 

and schematically is given in Figure 4. It is easy to see that the controller is of the type of a lead compensator 



Figure 4. Schematics of the LTI M-MRAC system. 


plus integral control, where the lead factor is determined by the parameter A. Therefore, A determines the 
phase margin (or the time delay margin) of the system, whereas 7 determines the performance. The ” sweet” 
spot given by the formulae ( 12 ) defines the relation between the design parameters that minimizes the input 
and output tracking errors. Increasing A from the ’’sweet” spot increases the robustness of the system, while 


decreasing it improves the performance. This is in contrast to conventional MRAC design (A = 0), where 
the lead compensator is absent, and the controller is pure integral controller. 

To get more analytic insight into the M-MRAC design, we consider a first order system with A rn = —a m , 
B = b > 0 and b m = a m . That is 


z(s) = — — u(s), x m = — — — r(s) 4 — [x(s) - x m (s)], u(s) = --[a;(s) - x m (s)\ , (31) 

s + a m s + a m s + a m s 


It is easy to see that the loop gain is 


L(s) = 


b 7 


s(s + a m + A) 

Therefore, the gain crossover frequency u> c can be computed from the equation 

67 


\ L (jUgc)\ = 


c \ / W? + (d m + A)"' 


= 1. 


which results in 


(32) 


(d m + A) 4 ^ ^2^,2 ( a ™ 


ujgc — 1 

The phase margin is computed from the equation 

PM = 7 r - arg (L(ju gc )), 

which results in 

/ 


(33) 


PM = tan 

2 


The corresponding time delay margin is 


\ 


l 

4 


b 2 y 2 


(dm + •^)‘ 1 


td = 


PM 


UJn 


(34) 


(35) 


It can be shown by direct computations that 


duj 


gc 


> 0 , 


dPM 


< 0 , 


dt d 

<9y 


< 0 


dj ’ dj 

That is, the phase margin and time delay margin decrease as 7 is increased for the MRAC and M-MRAC 
systems, which is consistent with nonlinear analysis. However, in the M-MRAC systems the desired robust- 
ness can be achieved with the proper choice of the parameter A, which is not available for MRAC systems. 
It can be shown by direct computations that 

du) gc dPM dt d „ 

d\ ’ dX > ’ dx > ■ 

Therefore, the time delay margin can be increased by increasing A. In fact, the desired time delay margin 
td can be maintained by setting A according to the equation 


cot 



(d m + A) 4 ^ ^2^,2 


N 


1 b 2 7 2 
4 (d m + A) 4 


1 

2 ' 


Figures 5 show the Bode plots for L(s) for a = 2, 6=1, 7 = 10000 in the A = 0 case, which corresponds 
to the MRAC design, and in A = 2^67 — a m case, which corresponds to the M-MRAC design. It can be 
observed that whereas the MRAC design has a phase margin of 1.15 degrees at gain crossover frequency of 
100 rad/sec, the M-MRAC design has a phase margin of 76.6 degrees at gain crossover frequency of 48.2 
rad/sec. 

Figure 6 shows the behavior of uj gc , PM and td with respect to A for 7 = 10000, when A is increased 
proportional to ^7. The decrease of u} gc , and increase of PM and td can be observed. 

We also provide the Bode plots of the closed loop system, to show that unlike the MRAC system, the 
M-MRAC does not exhibit any high gain behavior. As it can be seen from Figure 7, MRAC system’s gain 
has a high peak at the frequency y/b^y, but the gain of the M-MRAC system does not exhibit any peak. 



Bode Diagram 

Gm = Inf dB (at Inf rad/s) , Pm = 1.15 deg (at 100 rad/s) 


Bode Diagram 

Gm = Inf dB (at Inf rad/s) , Pm = 76.6 deg (at 48.2 rad/s) 




(a) MRAC (b) M-MRAC 

Figure 5. Bode diagrams with stability margins. 



Figure 6. Robustness properties with respect to A. 


Bode Diagram 

Gm = Inf dB (at Inf rad/s) , Pm = 1.62 deg (at 141 rad/s) 


Bode Diagram 

Gm = Inf dB (at Inf rad/s) , Pm = -1 80 deg (at 0 rad/s) 




Frequency (rad/s) Frequency (rad/s) 


(a) MRAC 


(b) M-MRAC 


Figure 7. Bode diagrams of the closed-loop systems. 


IV. Frequency Domain Analysis 

In this section we analyze the performance of the closed-loop M-MRAC system (31) in the presence of 
input and output disturbances. The system is displayed in Figure 8, where d is the input disturbance (load 



Figure 8. Closed-loop scalar M-MRAC system with input and output disturbances. 


disturbance) and n is the output disturbance (measurement noise). It is easy to see that the plant is 

b 


P(s) = 


s + a r j 


the controller is 


C(s) 


7 (s + am) 
s(s + a m + A) 


and the prefilter (reference model) is 


F(s) 


dm 

s + a m 


Therefore, the following relationships can be obtained from the systems block diagram [2] 


y 0) 

x(s) 

v(s) 

. «0 ) . 


r P(s)G(s)F(s) 

P(») 

1 



1 +P(s)C(s) 
P(s)C(s)F(s) 

1 +P(s)C(s) 
P(s) 

l+P(g)C(s) 

P(s)C(s) 


r(s) 

= 

l+P(s)C(s) 

C(s)F(s) 

l+P(s)C(s) 

1 

1 +P(s)C(s) 
C(s) 


d(s) 


1 +P(s)C(s) 
C(s)F(s) 

1 +P(s)C(s) 
P( s )C(s) 

l+P(s)C(s) 

C(s) 


. n ( s ) . 


L l+P(s)C(s) 

1 +P(«)C(s) 

1 +P(s)C(s) J 



which represent the transfer functions from the inputs r, d, n to the outputs x, y , u, v respectively. The 
six independent transfer functions involved in the above equation are called the Gang of Six in reference [2] . 
Their values are given in Table 3. 


Table 3. The Gang of Six for M-MRAC system. 


n <x\ n m 

n 1 n < ^ i a m 

~ ( s+0m )( s 2 + (0m+A)s+b7) 

— WrW ~ s 2 + ( 0m+A ) s+b7 

H d<i) n s(s+a m +A) 

n (cA n s(s+a m +A) 

^’yd\ b ) ^ xd{b) (s+am) ( s * + (am+x)s+b ^) 

<-x vd{b) s -2 + ( am+X )s+ b ' Y 

G xn (s) — G u d{s ) — s 2 + ( Qm+A ) s+f)7 

C ( „\ _ n t x\ _ 7(*+“m) 

^ un(K) — — s 2 + ( am+X ) s+bl 


First of all we notice that Theorem 2 gives the suboptimal value for A with respect to 7 according to 
equation (12). However, for this linear case we can select the optimal value for A as 


Ao = 2 Co-\At - a m , 

for the damping ratio Co = 1 / \/2 . In this case the phase margin (34) reduces to 


PM = — — tan 1 j * 


16Co 4 


(36) 


(37) 


which is independent of A and 7 and is equal to PM = 65.5, which is slightly lower than in the case of 
suboptimal A, as given in Figure 5. On the other hand, the gain crossover frequency depends on 7 and is 


expressed as 


“gc 


= V^7 


\ 



16C 0 4 


(38) 


For 7 = 10000, u gc = 64.4 rad/sec, which is slightly higher than in suboptimal case. The reason is that in 
the suboptimal the damping ratio is one. 

To check the control objective we apply the final value theorem to the transfer function from the distur- 
bance input d(s) to output x(s) 


T ry ( \ r s ( s + a m + A) 

= hm G x As) = Inn — — - 

s— >-0 s — >-0 (5 + tt m )(5 \CLm A)s ~\~ b 


= 0, 


where e ss is the steady state error due to a constant disturbance. If we select A = Ao, for large values of 7 
(fast adaptation) 


(s 2 + (a m + A )s -I- 67) 

and the closed-loop transfer function G xr {s) approximates the transfer function of the reference model F(s). 
Moreover, the transfer function from the command r to control u is approximately equal to 

®m7 ^ ®"rn 

( S 2 + (a m + A )s + 67) b 
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Figure 9. Bode diagrams for the linear M-MRAC system with the optimal value of A. 


This also can be observed from Figure 9, which shows the Bode magnitude plots for the transfer function 
Gy r (s), G ur (s) for different adaptation rates. Moreover, it can be seen that there is no high gain phenomena 
for large values of 7. 

To analyze the robustness of M-MRAC to input and output disturbances we plot the Bode magnitude 
for the remaining members of Gang of Six for different values of the adaptation rate, when the feedback gain 
A is computed according to equation (36). 

It can be observed from Figure 10 that all transfer functions are well behaved with the increase of the 
adaptation rate except for the noise sensitivity function G un (s). As Figure 10 (d) shows, the peak magnitude 
of the noise sensitivity function grows with 7 at higher frequencies. Therefore, the measurement noise will 
be amplified in the control channel. This phenomenon is similar to one reported in [9] for the conventional 
MRAC with cr-modification. However, MRAC with cr-modification exhibits a nonzero steady state error, 
thus failing to meet the control objective. 

Next, we fix 7 and vary feedback gain A. Figure 11 presents the Body magnitude plots for the members 
of the Gang of Six for 7 = 10000 and A increased from the optimal value Ao (36) by order of two. It can be 
observed that increasing A decreases the peak values of the noise sensitivity function, while the rest of the 
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Figure 10. Bode diagrams for the linear M-MRAC system for different adaptation rates. 


transfer functions do not experience significant changes. Hence, the robustness to measurement noise can be 
increased as desired by the proper choice of A. Meanwhile, the transfer functions from the command r(s) to 
the output y(s ) and control signal u(s ) remain in the acceptable range as it can be concluded from Figure 
12 . 


V. Robustness to Unmodeled Dynamics 

In this section we analyze the robustness of M-MRAC to unmodeled dynamics for an example of a first 
order system. The example is taken from reference [6] (p. 549). The system has a form 

y( s ) = ~~~ — (1 + A(s))u(s) , (39) 

s + a 

where a > 0 is a known constant, b > 0 is an unknown constant parameter, A(s) = — is the transfer 
function of the unmodeled dynamics, and /i > 0 is a small parameter. The objective is to track the output 
of the reference model 


n (s) , 


Vr{s) = 


(40) 
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Figure 11. Bode diagrams for the linear M-MRAC system for different values of A. 
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Figure 12. Bode diagrams for the linear M-MRAC system with different values of A. 


where r is bounded and piecewise uniformly continuous external command. The M-MRAC is designed for 
the reduced system (that is ignoring the unmodeled dynamics) in the form 

Vm(s) = -^—r(s) + —^—[y(s)~y m (s)} 

s + a s + a 

t K * ) = -7 [y(t) - ym(t)\r(t) 

u(t) = 0(t)r(t ) , 


(41) 


and is applied to the actual system, which results in the following error dynamics 


e(t) = —(a + A)e(i) + bz(t) — b6(t)r(t) — ar(t) 

z(t ) = ——z(t) — —9(t)r(t) 

Ht) = -7 e(t)r(t) . 

If r = constant , the system (42) is LTI with the matrix 


—a — A b — br 

0 -I 2r 

— yr 0 0 

the characteristic equation of which has the form 


s 


3 


a + A H — 


/ a + A 
V I 1 



b^r 2 

s + — — 

9 


= 0. 


( 42 ) 


Using the Routli-Hurwitz criterion on can show that the stability of the system depends on the sign of the 
third term of the first column of Routh table or equivalently on the sign of 


h = (a + A) 2 + ( — — bujr 2 ) (a + A) — 2 bjr 
9 


When A = 0 (conventional MR AC), one has h < 0 whenever 

2 1 a + a 2 n 

syr o> 

/i 2b + nab ' 

hence the characteristic equation has two roots in the open right-hand side s-plane, which coincides with 
the result derived in [6] in the case of a = 1 and b replaced with 2b — 1. That is, a fast adaptation in the 
conventional MRAC design leads to instability in the presence of the unmodeled dynamics A(s). Unlike 
MRAC design, M-MRAC design is robust to considered unmodeled dynamics, because we can make h > 0 
by setting 


A + a > — ^ + y - ^ 'bfijr 2 ^ + 2b r yr 2 . (43) 

In this case all roots of the characteristic equation are in open left-hand-side s-plane, hence the system (42) 
is exponentially stable, no mater how big y is. 

Figure 13 displays the M-MRAC performance from zero initial condition for a = 1, b = 2, r = 1 and 
y = 10000 (fast adaptation). A is set equal to the right hand side of the inequality (43), where b is replaced 
with a known upper bound bo = 5. It can be observed from Figure 13 (a) that in addition to asymptotic 
tracking, system’s both output and input signals have good transient response in the presence of unmodeled 
dynamics with n = 0.1. When the unmodeled dynamics are more significant (/i = 0.5), the M-MRAC still 
meets the control objective, but the transient behavior degrades, as it can be viewed in Figure 13 (b). 


VI. Conclusions 

We have presented robustness analysis of M-MRAC systems to time delays. The Razumiknin’s method 
was adopted to show that the time delay margin of the M-MRAC systems can be improved by increasing 
the error feedback gain, whereas the large adaptation rate improves the tracking performance. This tradeoff 
has been confirmed analytically via the frequency domain analysis in the case of linear adaptive systems. 
In the linear case the robustness and performance analysis have been also provided in the presence of load 
disturbances and measurement noise. The robustness of the M-MRAC design to unmodeled dynamics is 
analyzed for a first order system in comparison with the conventional MRAC system. 

The future research will address the robustness of M-MRAC design to measurement noise and unmodeled 
dynamics in general case of adaptive systems. 
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Figure 13. Bode diagrams for the linear M-MRAC system with the optimal value of A. 
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